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The radiative neutron capture cross section of 20r Pb has been measured at the CERN neutron 
time of flight installation n_TOF using the pulse height weighting technique in the resolved energy 
region. The measurement has been performed with an optimized setup of two CeD6 scintillation 
detectors, which allowed us to reduce scattered neutron backgrounds down to a negligible level. 
Resonance parameters and radiative kernels have been determined for 16 resonances by means of 
an R-matrix analysis in the neutron energy range from 3 keV to 320 keV. Good agreement with 
previous measurements was found at low neutron energies, whereas substantial discrepancies appear 
beyond 45 keV. With the present results, we obtain an s-process contribution of 77±8% to the solar 
abundance of 207 Pb. This corresponds to an r-process component of 23±8%, which is important for 
deriving the U/Th ages of metal poor halo stars. 

PACS numbers: 25.40.Lw,27.80.+w,97.10.Cv 



I. INTRODUCTION 



II. EXPERIMENT 



Since 207 Pb is one of the final products of the Th/U 
a-decay chains, its abundance provides a constraint for 
the Th/U abundances and their use as cosmochronome- 
ter 0]. The s abundance of 207 Pb has been shown to 
depend only weakly on details of stellar s-process models 
so that the r component can be reliably determined 
by subtraction from the solar value, N r = Nq — N s . 



Apart from its astrophysical importance, the neutron 
capture cross section of this isotope is also of relevance for 
the design of fast reactor systems. An eutectic mixture 
of Pb/Bi is presently considered as an appropriate spalla- 
tion target and coolant for accelerator driven systems . 
Because 22.1% of natural lead consists of 207 Pb, the rel- 
atively large neutron capture cross section of this isotope 
affects the neutron balance, and is therefore important 
for the design of this type of hybrid reactors. 



Given the large neutron scattering width in some of the 
207 Pb resonances, the neutron sensitivity of the detec- 
tor system becomes instrumental for the effective reduc- 
tion of backgrounds due to scattered neutrons. For this 
reason an optimized setup based on CqDs detectors has 
been employed by the n_TOF collaboration, well suited 
for the subsequent application of the pulse height weight- 
ing technique (PHWT), which was used in data analysis 
(Sec. IHIjl . With the adopted experimental setup angular 
distribution effects of the primary 7 radiation emitted 
for neutron capture events with orbital angular momen- 
tum I > are important. This effect, which is large for 
some resonances, could be properly treated in the data 
analysis as described in Section ITTT1 The resulting neu- 
tron capture cross section is presented in Sect ion llVl and 
the astrophysical implications of this measurement are 
summarized in Section Ivl 



The measurement was performed with an enriched 
sample in the form of a metal disk 20 mm in diame- 
ter and 2 mm in thickness, containing 92.40% of 207 Pb, 
5.48% of 208 Pb and 2.12% of 206 Pb. The sample was 
mounted in vacuum inside a sample changer made from 
carbon fiber together with a gold sample for absolute 
cross section normalization. In addition, a 208 Pb sample 
was used to study the background coming from in-beam 
7 rays, which are scattered by the sample. The long flight 
path of 185.2 m and the short proton pulse width of 6 
ns (RMS) are the important properties for achieving the 
high resolution in timc-of-flight (TOF) characteristic of 
the n_TOF spallation source [J. 

Neutron capture events were registered via the prompt 
capture 7-ray cascade by a set of two CeD 6 detectors, 
which were optimized with respect to neutron sensitiv- 
ity p|. The detectors were placed at ^125° with respect 
to the incident neutron beam in order to minimize angu- 
lar distribution effects of the emitted capture radiation. 
A schematic view of the experimental setup can be seen 
in Fig. 2 of Ref. In this configuration, the in-beam 
7-ray background was also substantially reduced. 

The neutron flux was monitored by means of a 
200 /ig/cm 2 thick 6 Li-foil mounted 3 m upstream of the 
207 Pb sample. Particles from 6 Li(n, cv) 3 H reactions are 
registered with four silicon detectors surrounding the 6 Li- 
foil outside of the beam |7| . 

The saturated resonance technique Q using the 4.9 eV 
197 Au resonance was applied for absolute calibration of 
the 207 Pb capture yield. Calibration of the yield in this 
way requires the precise knowledge of the neutron in- 
tensity versus the neutron energy. This has been deter- 
mined with an uncertainty of 2% [tj by two independent 
measurements performed with the 6 Li monitor described 
before and with a calibrated fission chamber ^1 ■ 



III. CAPTURE DATA ANALYSIS 



'Corresponding author, address: Apdo. Correos 22085, 
46071 Valencia (Spain). Tel.:+34963543499, e-mail ce- 

sar.domingo.pardo@cern.ch. 



The first aspect of the data analysis obviously concerns 
the determination of the weighting function (WF). Based 
on previous experience, the WF for the measured sam- 
ples of 197 Au (used for normalization) and of 207 Pb was 
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calculated via the Monte Carlo technique. The procedure 
followed has been described in detail in Refs. pdlll2|. 

The capture 7-ray spectra of both samples were calcu- 
lated with a Monte Carlo code in order to estimate the 
corresponding uncertainty of the WF, which turned out 
to be less than 0.5%. It has been experimentally demon- 
strated ^3 that the combination of WFs obtained by the 
Monte Carlo method with the saturated resonance tech- 
nique yields an overall systematic uncertainty of better 
than 3%. However, this level of accuracy can be only 
achieved if all sources of systematic uncertainty are prop- 
erly taken into account. This refers mainly to (i) the 
determination of the neutron flux, (ii) the treatment of 
the background components, and (iii) the effect of the 
electronic threshold used for the signals from the CgD 6 
detectors. In the particular case of 207 Pb, one has to 
consider also a pronounced angular distribution effect of 
the capture 7 rays due to the low multiplicity in the de- 
excitation pattern of 208 Pb. This effect will be considered 
separately since it depends strongly on the resonance spin 
and parity (J 17 ). 

The various sources of systematic uncertainty are listed 
in Tabled In the following, the treatment of these effects 
is described in detail. 

TABLE I: Systematic effects and related uncertainties. 
Effect Uncertainty (in %) 



Weighting function, saturated resonance 

technique and electronic threshold 
Background determination 
Energy-dependence of neutron flux 
Neutron sensitivity 
Angular distribution: for J n — 1 + 
for r=2+ 



<2 

<0.5 (1.5) a 
2 

<0.3 
0.6 (8) 6 
0.3 



a Broad s-wave resonances at E a = 41 and 256 keV (see also 
Table EJ. b Resonances with unknown angular distribution. 



for determination of the background from in-beam 7 rays. 
The respective systematic uncertainties were less than 
1.5% for the broad resonances at 41 and 256 keV, and 
below 0.5% for all the others. 

A different type of background might arise in the mea- 
surement of resonances that show a dominant scatter- 
ing channel like the s- and p-wave resonances at 41 keV 
and 128 keV (Table EJ, where r„/T 7 w 300. Thanks 
to the optimized detection setup and the small amount 
of material around the sample and the detectors, this 
background turned out to be negligible in the present 
measurement. 



B. Electronic threshold 

The low energy cut-off in the pulse height spectra due 
to the electronic threshold («340 keV in this measure- 
ment) had a non-negligible influence on the yield mea- 
sured with the 7-ray detectors. The effect can be cor- 
rected by modeling the capture cascades with a Monte 
Carlo simulation of the complete pulse height spectra 
recorded in the n_TOF detection system as described in 
Refs. HE1 El- 
Following neutron captures on 207 Pb, the de-excitation 
pattern of the 208 Pb resonances is rather simple. It ba- 
sically consists of a single 7-transition of 7.37 MeV for 
J = 1 and of a two-step cascade in case of J = 2 reso- 
nances Therefore, the Monte Carlo simulations 
of the capture spectra and the estimate of the threshold 
effect are correspondingly straightforward and reliable, 
yielding correction factors of 5±1% and 4±1% for J = 1 
and J = 2 resonances, respectively. 



C. Angular distribution effects 



A. Background 

In the present measurement the background in the re- 
solved resonance region (RRR) is dominated by delayed 
7-rays accompanying the neutron beam and scattered in 
the sample. These 7-rays arise essentially from neutron 
capture in the water moderator of the n_TOF spallation 
target and exhibit a smooth dependence on neutron time 
of flight. For samples of equal atomic number Z and 
similar thickness the background level is also very sim- 
ilar. Since all observed resonances in 207 Pb were well 
isolated, the background level was best determined by 
choosing a relatively wide neutron energy window around 
each resonance and defining the background level by a 
constant term. The neutron energy dependence of the 
background level fitted for each resonance was afterwards 
cross checked with a measurement of a 208 Pb sample. 
The 208 Pb sample itself showed very few resonances in 
the entire energy region and was, therefore, well suited 



The low multiplicity (m = 1 for J = 1 and m = 1- 
2 for J = 2 resonances) of the 207 Pb capture cascades 
has the experimental disadvantage that it causes strong 
angular distribution effects in the measured pulse height 
spectra. Indeed, captures with orbital angular momenta 
I > lead to aligned states in the compound nucleus, 
perpendicular to the incident neutron beam. This causes 
anisotropy in the angular distribution of the prompt 7 
rays, 

W{0) = A k P k {cos0) = 1 + A 2 P 2 (cos6) + 

k 

+A A P 4 (cos9) + A 6 P 6 (cos9) 7 (1) 

where Pk(cos9) are the Legendre polynomials of order 
k and A k are coefficients, which depend on the initial 
(J) and final (</') spin values, on the multipolarities (L) 
of the transition, and on the degree of alignment. With 
ideal 7 detectors of negligible volume, this effect would be 
minimized by setting both detectors at 125°. However, 
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the CeD6 detectors (~1 1 volume) cover a substantial 
solid angle. Hence, 7-rays are registered in a relatively 
broad angular range around 125°. 

For resonances with spin J n = 1 + there is a channel 
spin admixture (s = 0,1), which contributes to an in- 
complete alignment with generally unknown proportions 
of the two channels with s = and s = 1. The A/, co- 
efficients in Eq. ^ could be determined for the two 1 + 
resonances at 30.5 keV and 37.7 keV (Table[nJ by means 
of measured angular distributions |l5| and using the for- 



mulae in Ref. [16J . With this information the Monte Carlo 
simulations of the experimental setup yielded correction 

factors of f^Q keV = 0.965(3) and f^^eV — 1.037(6) as 
described in detail in Ref. For the other three 1 + 
resonances at 90 keV, 128 keV and 130 keV, where angu- 
lar distributions are unknown, the corrections could not 
be determined. Instead, a systematic uncertainty of 8% 
was adopted in these cases. This value was estimated by 
means of Monte Carlo simulations of our experimental 
setup using different angular distributions, ranging from 
alignment zero up to total alignment. 

For resonances with J n =2 + the de-excitation occurs 
predominantly by emission of a two-step cascade [l3j . 
The angular distribution of the first 7-ray can be cal- 
culated using the formulae in Ref. [16|. The angular 
distribution of the second 7-ray is influenced by the re- 
alignment introduced by the previous transition and it 
could be also calculated by modifying the formulae in 
Ref. (lt| . The final yield correction factor for 2 + reso- 
nances (considering also a 10% branching to the ground 
state) results in a value f-" 1 = 1.015(3). 




3.04 3.06 3.08 3.1 

E n (keV) 




10.3 
E n (keV) 



FIG. 1: (Color online) R-matrix fit for two of the resonances 
measured in this work at 3 keV and 10 keV. 



D. R-matrix fits 

The capture yield corrected for the effects described 
above, 

f x f xY ex P = B + Y(E ,r i ,r n ), (2) 

was analyzed by means of the R-matrix analysis code 
SAMMY • -B is a constant term describing the back- 
ground in the region of each resonance. In cases, where 
the neutron width T n is well known from transmission 
measurements and considerably larger than the capture 
width r~, r„ was kept fixed in our analysis. In other 
cases and where the number of counts in the resonance 
was sufficiently large, we preferred to vary r„ and T 7 in 
order to better describe the corresponding capture area 
or radiative kernel. Fig. ^ shows the measured capture 
yield for the two first resonances in 207 Pb. The contin- 
uous line corresponds to an R-matrix fit performed with 
the SAMMY code. 

IV. RESULTS AND UNCERTAINTIES 

The parameters and radiative kernels of the analyzed 
207 Pb resonances are summarized in Table |H] Orbital 



angular momenta I and resonance spin J were taken from 
Ref. [3. 

According to the discussion in the previous section, the 
radiative kernels can be given with an overall systematic 
uncertainty of «3%, except for the resonances at 90 keV, 
128 keV and 130 keV with spin J = 1, where a systematic 
uncertainty of 8% had to be adopted (see Table P). 

The radiative kernels (last column in Table [H| are 
compared in Fig . [21 with previous measurements made 
at ORNL[H[l||. 

For the first four resonances our results are in good 
agreement, whereas serious discrepancies appear at 
higher energy, mainly beyond 30 keV. It is remarkable 
that the radiative kernels of seven out of the nine J = 1 
resonances reported here are systematically smaller than 
those of Ref. [lj] , whereas the cross sections for J = 2 res- 
onances are in agreement or higher. Such discrepancies 
can not be explained only in terms of the experimental 
aspects discussed in the previous section, but must be 
due to the hardness of the capture 7-ray spectra. Since 
,/ = 1 resonances exhibit substantially harder spectra, 
the discrepancy suggests a problem with the WF used in 
the previous work. A further hint in this direction is ob- 
tained from the fact that a thin sample (0.5 mm in thick- 
ness) was used in Ref. [l3l | for the neutron energy range 
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TABLE II: Resonance parameters and radiative kernels from 
the analysis of the 207 Pb(n,7) data measured at n_TOF a . 



E 


/ 


./ 


r„ 


r 7 


,gr 7 r„/r 


(eV) 






(mev) 


(mev) 


(mev) 


3064.700(3) 


1 


2 


111.0(8) 


145.0(9) 


78.6(9) 


10190.80(4) 


1 


2 


656(50) 


145.2(12) 


149(14) 


16172.80(10) 


1 


2 


1395(126) 


275(3) 


287(30) 


29396.1 


1 


2 


16000 


189(7) 


234(9) 


30485.9(5) 


1 


1 


608(45) 


592(50) 


225(30) 


37751(3) 


1 


1 


50xl0 3 


843(40) 


620(30) 


41149(46) 





1 


1.220xl0 6 


3970(160) 


2970(120) 


48410(2) 


1 


2 


1000 


230(20) 


235(20) 


82990(12) 


1 


2 


29x 10 


360(30) 


A A A { f\C\\ 

444(30) 


90228(24) 


1 


1 


272 xlO 3 


1615(100) 


1200(80) 


127900 


1 


1 


613xl0 3 


1939(150) 


1449(120) 


130230 


1 


1 


87xl0 3 


900(80) 


675(60) 


181510(6) 





1 


57.3xl0 3 


14709(500) 


8780(300) 


254440 


2 


3 


lllxlO 3 


1219(90) 


2110(150) 


256430 





1 


1.66xl0 6 


12740(380) 


9482(280) 


317000 





1 


850 xlO 3 


10967(480) 


8120(350) 



Orbital angular momenta I and resonance spins J are from 
Ref. 0. 
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FIG. 2: (Color online) Comparison of the radiative cap- 
ture kernels for 207 Pb resonances measured at n_TOF and 
ORNL [T3I H4|. In the bottom panel, the ratio between the 
two data sets is shown. 



uncertainties, both cross sections are generally in good 
agreement. However, the new MACSs are clearly higher 
at energies below kT =15 keV. In this region, the un- 
certainties could be reduced from 13% down to 5%, an 
improvement by more than a factor of two. 
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FIG. 3: (Color online) Maxwellian averaged cross sections for 
a range of thermal energies compared to previous values [l9| . 

In the calculation of the MACS our results were 
complemented with some additional resonances from 
Ref. [23, which could not be observed at n_TOF due 
to the in-beam 7-ray background. However, the contri- 
butions of these resonances are rather small, i.e. only 2% 
and 7% of the MACS at kT = 5 and 25 keV, respectively. 



TABLE III: Maxwellian averaged cross sections of 207 Pb. 



kT 
(keV) 



Ref. [19] 



MACS (mb) 



This work 



10 

20 
25 
30 



13.9 (19) 

12.3(16) 
11.5 (16) 
10.7 (14) 

9.7(13) 



16.2(8) 
14.1(8) 
13.7(9) 
12.0(7) 
10.9(6) 
9.9(6) 



below 45 keV, whereas a much thicker sample (16 mm) 
was employed above that energy. The better agreement 
with the thin-sample measurement su gges ts that the WF 
calculated for the analysis of Ref. [l3|, [LJ| was only valid 
for the thin sample case, reflecting the strong influence 
of the s amp le thickness on the shape of the WF as noted 
in Ref. [TJ. 

V. IMPLICATIONS FOR THE s- ABUNDANCES 
IN THE PB-BI REGION 

The Maxwellian averaged cross section (MACS) ob- 
tained from the present results are compared in Fig. |3| 
with the values reported in Ref. [nj . Within the quoted 



According to the Galactic chemical evolution (GCE) 
model described in Refs. 0, |23, the s-process abun- 
dance of 207 Pb is essentially produced in low mass asymp- 
totic giant branch (agb) stars. In these AGB stars energy 
is produced in a thin layer surrounding the inert C/O 
core. Subsequently, extended and quiescent H-shell burn- 
ing episodes are followed by much shorter He shell flashes 
[23. In this model about 95% of the neutron exposure is 
due to the 13 C(a,n) 16 reaction, which operates during 
the interpulse phase between He shell flashes at tempera- 
tures around ~ 10 s K, corresponding to a thermal energy 
of kT w 8 keV. At this stellar temperature the present 
MACS is about 13% higher than the value reported in 
Ref. In order to estimate the effect of the higher 

cross section on the calculated s-process abundances, a 
model calculation was made for thermally pulsing AGB 
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stars with M = 3M©, M = 1.5M© and a combination of 
metallicities: [Fe/H] = —0.3 (most characteristic of the 
main s component [24|). and [Fe/H] = —1 (representa- 
tive of the strong s component 0>|2j|). In this way, the 
s-process fraction of 207 Pb could be determined as N s 
= 77(8)%, instead of 82(18)% 0. In the present deter- 
mination of N s for 207 Pb, the contributions obtained for 
the main and strong components were 60% and 17% of 
the solar 207 Pb abundance, respectively. In the evalua- 
tion of the final uncertainty of the s-process abundance 
of 207 Pb, the estimated contribution due to the uncer- 
tainty in its cross section is now small, of about 3%. The 
contribution due to the uncertainty in the neutron cap- 
ture cross sections of 204 Pb, 205 Pb, and 206 Pb [HI has 
been estimated as less of 2%. The uncertainties of the s- 
process model are estimated to be ±3% for the main and 
±10% for the strong component, resulting a contribution 
of less than 4% to the total uncertainty in the s-process 
abundance of 207 Pb. Finally, the main source of uncer- 
tainty in the determination of N s for 207 Pb is due to the 
uncertainty in the solar abundance of lead, which is of 

8-10% jH mil. 

Since the s abundance of this isotope is robust with 
respect to uncertainties in the s-process model, the r 
component can now be constrained as N r — 23(8)%. 
Within the quoted uncertainty of 35%, this result is 
compatible with earlier r-process calculations 1| based 
on the ETFSI-Q nuclear mass model, which find an r- 
process fraction for 207 Pb of 18.4% with an uncertainty 
of 10-20%, and with r abundance calculations reported 
more recently |2g, which yield values between 15.1% and 
16.4%. 



resulting in an r fraction of 23(8)%. 
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VI. CONCLUSIONS 



The neutron capture cross section of 207 Pb has been 
measured at the high resolution neutron time-of-flight 
facility n_TOF. 

The main sources of systematic uncertainty affecting 
this measurement have been treated in detail. From 
the experimental point of view background due to scat- 
tered neutrons could be eliminated by means of an op- 
timized detection setup, and angular distribution effects 
were minimized by setting the detectors at 125° with 
respect to the neutron beam. In data analysis, the re- 
maining systematic effects have been carefully corrected, 
resulting in a 3% accuracy for the radiative kernels of the 
capture resonances, except for three cases, where missing 
angular distribution data led to an uncertainty of 8%. 

We report resonance parameters and capture areas for 
16 resonances in the neutron energy interval from 3 keV 
to 320 keV. At low neutron energies our results are in 
good agreement with previous data, but reveal significant 
discrepancies for resonances above 40 keV. Maxwellian 
averaged cross sections for 207 Pb were determined with 
an accuracy of ±5%. With this information the s-process 
component of solar 207 Pb was determined to be 77(8)% 
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